Single-ion anisotropy in Haldane chains and form factor of the 0(3) nonlinear sigma model 
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We consider spin-1 Haldane chains with single-ion anisotropy, which exists in known Haldane chain materi- 
als. We develop a perturbation theory in terms of anisotropy, where the magnon-magnon interaction is important 
even in the low temperature limit. The exact two-particle form factor in the 0(3) nonlinear sigma model leads to 
quantitative predictions on several dynamical properties, including the dynamical structure factor and electron 
spin resonance frequency shift. These agree very well with numerical results, and with experimental data on the 
Haldane chain material Ni(C5Hi4N2)2N3(PF(i). 
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One-dimensional quantum spin systems are an ideal subject 
to test sophisticated theoretical concepts against experimental 
reality.— One of the best examples is the Haldane gap problem. 
Haldane predicted in 1983 (Ref.i) that the standard Heisen- 
berg antiferromagnetic (HAF) chain 'H - J J2j ' ■^j+i has a 
non-zero excitation gap and exponentially decaying spin-spin 
correlation function for an integer spin quantum number S . 
It has been long known that the HAF chain with S - 1/2 
is exactly solvable by a Bethe ansatz, and that it has gapless 
excitations and the power-law spin-spin correlation function. 
While the same model cannot be solved exactly for S > I, 
Haldane's prediction was rather unexpected and surprising at 
the time. 

Haldane's argument was based on the mapping of the HAF 
chain to the 0(3) nonlinear sigma model (NLSM), which is a 
field theory defined by the action 
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where g = 2/S is coupling constant, v is spin- wave velocity, 
- InS and Q = (l/47r) J dtdx n ■ dtn x dxti is an integer- 
valued topological charge. The field n{x) is related to the 
spin Sj via Sj ~ (-1)-' ylS{S H- 1) «(x) -H L{x), where L{x) - 
n X d,n/g. The field n has a constraint = 1. For a half- 
integer S , the topological term iOQ should be kept. However, 
for an integer S , the topological term iOQ = 2m x (integer) 
is irrelevant and it suffices to drop iOQ in eq. ([T]). The 0(3) 
NLSM without the topological term is a massive field theory, 
which implies that the integer S HAF chain (Haldane chain) 
has a non-zero gap and a finite correlation length. The Hal- 
dane's conjecture is now confirmed by a large body of the- 
oretical, numerical, and experimental studies.- Moreover, the 
0(3) NLSM is also useful in describing integer S HAF chains. 

There are various complications in real materials. A Hal- 
dane chain material generally has a single-ion anisotropy 
(SIA): -H' = Y.jlD{S'f + E\{S])^ - (S-j)^]]. This interaction 
is important, for example, for electron spin resonance (ESR) 
measurements. ESR is a useful experimental probe which 
can detect even very small anisotropics. In other words, the 
anisotropic interaction is the key to understanding a rich store 
of ESR experimental data. However, the theory of ESR is not 



sufficiently developed for many systems, including Haldane 
chains, leaving many experimental data not being understood. 
In order to fully exploit the potential of ESR, accurate formu- 
lation of the SIA in Haldane chains is required. 

The SIA can be treated as a perturbation since it is usually 
small compared to the isotropic exchange interaction /. In the 
0(3) NLSM language, the perturbation is written as 



"K' = 5(5 4- 1) J dx[D(n~f - E{(n'f - (n>f 
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which spoils the integrability of the 0(3) NLSM. Several 
simple calculations have been done based on the Landau- 
Ginzburg (LG) model .iii. When the elementary excited par- 
ticles (magnons) are dilute, the interaction between magnons 
may be ignored. If this is the case, the system is effectively 
described by a much simpler theory of free massive magnons 
(the LG model).— However, description by the LG model is 
not accurate and, furthermore, it is phenomenological.— Even 
in the low-energy limit, where the free magnon approximation 
is supposed to be exact, it is not the case with respect to the 
evaluation of Eq. (|2]i. This is because the perturbation (|2| cre- 
ates and annihilates two magnons at the same point; in such a 
situation, interaction among the magnons is indeed important 
even when the average density of magnons in the entire sys- 
tem is infinitesimal. Therefore, correct handling of the SIA in 
the 0(3) NLSM framework requires a proper inclusion of the 
magnon interaction. 

In this paper, we present such a formulation, utilizing the 
integrability of the 0(3) NLSM. The effects of interaction are 
encoded in the form factors of operators. The form factors 
in integrable field theories can be determined by the consis- 
tency with the exact 5 -matrices and several additional ax- 
ioms i^"— Form factor expansion (FEE) is particularly power- 
ful in massive field theories such as the 0(3) NLSM, because 
the higher-order contributions survive only above the higher 
energy thresholds.'- The leading contribution to the FEE of 
Eq. (|2]i is given by the two-particle form factor. The FEE 
shows an excellent agreement with the correlation function of 
(S ')^ numerically obtained in the 5=1 HAF chain, demon- 
strating the importance of the interaction. At the same time, 
the renormalization factor for the SIA (|2]i is determined by the 



fitting of the numerical data. Furthermore, we discuss two ap- 
pHcations to physical problems of interest: the split of triplet 
magnons in the dynamical structure factor and the ESR shift in 
the 5 = 1 HAF chain with SIA. We find very good agreement 
with numerical results in both applications, and with experi- 
mental data on the ESR shift, without introducing any extra 
fitting parameter 

A single magnon excitation can be parametrized by the 
rapidity 0, so that its energy and wavenumber are given 
respectively as Aocoshfl and (Ao/v)sinh0, where Aq = 
0.417 is the Haldane gap. Because of interactions among 
magnons, the S matrix of 0(3) NLSM has a complicated 
structureJi The one-particle form factor of an operator O 
is defined as a matrix element which connects the ground 
state |0) to a one-particle state \Oi,a\) (a\ - 1,2,3), 
namely Foi6\,ai) = {<d\0\9\,ai). And the «-particle 
form factor is defined as Fo{0\,a\\92,a2\- ■ ■ ',d„,a„) = 
{Q\0\9i,a\\62,a2\- ■ ■ ;9n,aii), where this «-particle state 
is normalized as {9'y,a\;---;9'^,a'J^6\,a\\---;9„,a„) - 

The FFE of the fundamental field n", which corresponds to 
(a staggered part of) the spin operator 5'', has often been stud- 
ied. The leading contribution to the FFE is the one-particle 
form factor F„ci{9\,a\). Because n" is odd under the trans- 
formation n — > -n, the next order contribution comes from 
the three-particle form factor, which gives small corrections to 
the spin-spin correlation function. '^ '^ On the other hand, the 
composite operator {S ")^, which is of our central interest, has 
been less studied. Since it is proportional to («")^ and even un- 
der the reversal n -n, the leading contribution to the FFE 
comes from the two-particle form factor F(jjay2{6i,a\; 62,02) ■ 
We note that the exact two-particle form factor of the anti- 
symmetric field L(x) in the 0(3) NLSM has been applied to 
describe the uniform part of the spin-spin correlation function 
of HAF chainsJ^— Including the renormalization factors for 
spin operators, which are undetermined at this point, we have 

Fs40i,ai)^ ^/Z6a,a„ (3) 

F(s»)2(0i,fli;02,fl2) = -iZ2 6a„a,_(36a,a, ' 1)«A2(01 - O2). (4) 

The two-particle form factor (HJi receives contributions from 
higher-order terms in the FFE ofS", and cannot be determined 
by Eq. (|3]l alone. Thus Z2 is a parameter independent of Z. 

We have the constraint X]fl=i 2 3^'^°)^ = 2 on the compos- 
ite operator. From this constraint and the 0(3) symmetry, it 
follows that E«=i,2,3^(S'O=(fi>3;02,3) = (Ol^i, 3; 02, 3> = 0, 
which is satisfied by (|4]l. Integral representation of i/'2(0) is 
given in Ref. |H, for O(A^) NLSM with a general integer A^. 
For = 3, it reads 
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Figure 1. (color online): Numerically calculated spin-spin 
correlation (- l)''(0|5'(r)5'(0)|0> (circles) and the correlation 
<0|(5~(f-))2(5'(0))2|0> - 4/9 (triangles) are compared with FFEs © 
with Z = 1.26 (solid curve) and the connected part of Q with 
Z2 = 0.24 (dashed curve). The free magnon approximation (dotted 
curve) cannot fit the correlation function of (5~)^. 



the FFE for {S")^ and further to determine Z2, we computed 
the equal-time correlation function {Q\{S\r))\S^{0)f\0) 
by the infinite time-evolving block decimation (iTEBD) 
method<2!i as shown in Fig. [T] 

FFE is derived by inserting the identity 1 = X^/'ilo^"' 
where the P„'s are the projection operators to the n-particle 
subspace of the Fock space, defined by Pq - |0)(0| and 



-\6i,au--- ■,0„,a„}{9i,ai; - ■ ■ ;6'„,fl„| 



for n > 1 . In the leading nonvanishing order, we find 
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Z = 1 .26 was given in RefiJ^ by comparing numerically ob- 
tained spin-spin correlation function with the LG model. Con- 
cerning the spin-spin correlation function, the LG model is 
equivalent to the lowest-order FFE ( |6]i; our iTEBD calcula- 
tion also reproduces the result of Ref. [l6l On the other hand, 
to the best of our knowledge, Z2 has not been determined pre- 
viously. 

As shown in Fig. [1] the lowest order of FFE (|7]i shows an 
excellent agreement with the numerical data; the fit also de- 
termines 



1/^2(6) = sinh - exp 
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^cosh[(7r + i9)Lij] - 1 



/o vj sinh(7rcLi) 
This integral can be analytically carried out to give 

«A2(0) = ^(0-^Otanh^. 



(5) 



Determination of the renormalization factors Z and Z2 re- 
quires numerical calculations. In order to test the validity of 



Z2 = 0.24. 



(8) 



Since we used the known values of the Haldane gap Aq = 

0. 417 and the spin-wave velocity v = 2.49/ (Ref.H]]) for S - 

1, the renormalization factor Z2 is the only fitting parameter. 
In contrast to the FFE (|7]i, the LG model, which ignores 

interaction among magnons, shows discrepancy with the nu- 
merical data, as also shown in Fig. [1] To illustrate the effect 
of the interaction, let us discuss the asymptotic long-distance 
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Figure 2. (Color online): Numerically determined excitation gaps 
Av (circles) and A- (triangles) are plotted for -0.4 < DIJ < 0.6 
and E = Q. Deviation ofthe numerical data from the first order 
FFPT (solid and dashed lines) is attributed to higher-order perturba- 
tions. Inset: The ratio <S~-(;r, Aj)/<S"(;r, Aj) obtained by the Lanc- 
zos method (symbols) and l ll4b (solid curve) are compared. The 
extrapolation to L = oo is done by fitting the finite-size data for 
L = 12, 14, 16, 18 and 20 with a polynomial of 1/L. 



behavior of Eqs. ^ and O. When r — > H-oo, only the be- 
havior of (/'2((?) at (? ~ is relevant in (|7]i. Here we can ex- 
pand © and © as (-l)'X0|5=(r)5~(0)|0) oc e-''^/ VS^iT/l and 
<0|(5-'(r))"(5-'(0))'|0> - 4/9 oc e-''^'- /(4nr/^2)- In a relativis- 
tic field theory, the inverse correlation length is equivalent to 
the lowest excitation energy created by the operator; in fact 
^ - v/Aq. Furthermore, in the LG modelr. ^2 - should 
hold. This is because the composite field (n")^ creates two par- 
ticles, and 0(3) NLSM does not contain any bound states.— 
Thus the excitation energy for the two-particle creation would 
be twice the magnon mass (2Ao), implying ^2 - ^/2. How- 
ever, the actual numerical data are inconsistent with this rela- 
tion: ^2 - 2.75 < ^/2 = 3.01. This discrepancy is attributed 
to the interaction between magnons. Since (n")- creates two 
magnons at the same point, the actual excitation energy is 
larger than 2Ao, resulting in ^2 < ^/2. 

With the full determination of the two-particle form factor 
(|4|l, we turn to discussion of dynamical structure factor (DSF) 
at r = in Haldane chains with a SIA. The peaks in the DSF 
reflect the energy of the magnon at a given momentum. Triply 
degenerate magnon dispersions in the isotropic chain are split 
due to the SIA. We determine the first-order perturbation to 
the masses. A*," = A„ - Aq, in the form-factor perturbation 
theory (FFPT)2^: 



{e,aYH'\Q,a) 
{6,a\Q,a) 



(9) 



In fact, both the numerator and the denominator are propor- 
tional to 6{0), and Eq. (|9]l should be understood as the ratio of 
the coefficients of 6{0). Furthermore, the numerator equals to 
F'j^>{Q, a; O-ni, a) because of the crossing symmetry There- 
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Figure 3. (Color online): Magnetic field dependence of ESR shift 
Yd{T,H) for r = O.iy (circles) and T = 0.27 (triangles). The solid 
curve is l |18i l. which is exact in // — » 0. The dashed and dotted curves 
are l |19t at T = 0.17 and T - 0.27, respectively. 



fore, (|9| reads 
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The leading contribution to the T - Q DSF S""{7t, lo) cor- 
responds to the creation of a single magnon. Therefore we 
find 



(13) 



which has the identical form to the DSF of a system of free 
particles. This is natural because the population of magnons 
approaches zero in the T — > limit, and thus the interactions 
are negligible. Nevertheless, we emphasize that the change of 
the masses as A^ (fT0i)- (fT2l) due to the SIA is affected by the 
magnon-magnon interaction. Equation ( fTJi l implies that the 
magnon masses A^, can be identified with the peak frequency 
of DSF at the antiferromagnetic wavevector q - 71. In Fig.|2] 
we compare the magnon masses A„ extracted from the T - Q 
DSF peak obtained numerically by the Lanczos method^^ for 
various values of D (while setting E - 0). For small D, the 
numerical data agree very well with the FFPT (fT0b-(fT2b. 

The form of the T - Q DSF ( fT3] ) leads to another prediction: 
The ratio of the DSF intensities should obey 



/ doj S'~{n, A,) 
/ doj S-'^in, A,) 



A, 
A/ 



(14) 



This is also confirmed by the Lanczos data as shown in the 
inset of Fig.|2] 

Let us extend our discussion to the system under a finite 
magnetic field. Now our Hamiltonian 'H - 'Hq + "Hz + "K' 
consists of three terms. Wo is the SU(2) symmetric exchange 
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Figure 4. (Color online) Comparison of the resonance frequency 
w,. = ge/Jfi// + 5ii> by QMC (circles) with experimental dat^ (tri- 
angles). We performed QMC calculations with L = 30 sites. We 
used D = 0.25 J and H \\ c {& = 9 = 0). The solid curve is obtained 
from ( 1191 ) and the dashed line represents the paramagnetic resonance 



interaction, 'Kz = -gefisH ■ S - -geHsH ■ J2i^i ^^e Zee- 
man interaction, and 'H' is the SIA, which is assumed to be 
small, ge is Lande g factor of electrons and yUs is the Bohr 
magneton. We set g^jJ^B = 1 unless otherwise stated. ESR 
is a very powerful tool to study the effects of anisotropics on 
spin dynamics. One of the fundamental quantities in ESR is 
the resonance frequency shift (ESR shift). The ESR shift is 
generally given, in the first order of the anisotropy 'H', as^"— 



{[['H',S%S-])o 



(15) 



(■ ■ ■ )() denotes the average with respect to the unperturbed 
Hamiltonian = "Ko + 'Hz- For the SIA, dllll reads 

Sto = /(©, O) Yd{T, H), where /(©, <5) = D(l - 3 cos- 0) - 
3E sin" © cos 2<1) and 



ZjOiSf -2)0 



(16) 



(©, <I>) is the polar coordinate of the magnetic field axis. 

To apply the results of the FFPT, first we consider the limit 
r, // «c Aq. Here we could project the numerator to one- 
magnon subspace, ignoring the multi-magnon contributions. 
The projection operator is Pi = / ^ ^a=o ± 1^' '^)^^' Note 
that we introduce a different set of indices a - Q, + represent- 
ing magnons with dispersion Ea(ff) = Aocosh^ - aH. The 
projection leads to 



Pi ^[3(5^)2 -2] Pi 



do 3Z2V 



An 2Ao cosh 6 
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Its thermal expectation value can be given in terms of the 
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Figure|3]shows the magnetic field dependence of Yd{T, H), 
comparing ( fTSl l from the FFPT with the numerical results ob- 
tained by ( fT6] l with quantum Monte Carlo (QMC) method in 
ALPS software.i^ 

Although the agreement is good at low temperature T = 
O.iy and at low magnetic fields H <si Aq, the discrepancy is 
evident for// > Aq. This is rather natural, because the magnon 
population increases as H is increased, invalidating the dilute 
limit approximation made in the derivation of Eq. ( fTSl ). In 
particular, T = 0, // = Aq is a quantum critical point which 
separates the low field gapped phase and the high field TLL 
phase, where magnons are condensed. Although it is difficult 
to handle the case with nondilute magnons, a reasonable im- 
provement would be incorporating magnon-magnon repulsion 
through the Pauli exclusion principle by utilizing the Fermi- 
Dirac distribution function fa{k) = [e'^'M/^ + i]-' instead of 
the classical one, in Eq. ( fTST l. This is demonstrated by the 
fact that the z - 2 free-fermion theory well describes the low- 
energy behavior near the quantum critical point H = Ao3^ 
The magnetization is = m(T,H) = / 5^ [f+{k) - f-(k)] 
and Yd{T, H) is 



Yd(T,H) 



3Z2 



dk 



2m{T, H) J In 2coo{k) 



[2Mk)-Mk)-Mk)]. 

(19) 



This reduces to Eq. (fTSl ) in the limit H,T — » 0. We em- 
phasize that there is no free parameter in our theory since 
the renormalization factor Zo in the overall coefficient of (fT9] l 
has been already determined in dH). As shown in Fig. [3] the 
free-fermion approximation (fT9] l explains the extremum of 
the ESR shift observed numerically around the critical field 
H = A,,. 

Figure |4] shows the ESR shift observed experimentally in 
Ni(C5Hi4N2)2N3(PF6),^- which possesses the SIA, and the 
corresponding numerical result by the QMC method. Our 
FFPT ( fT9] l successfully accounts for the experimental and nu- 
merical results, including the gradual approach to the param- 
agnetic resonance line co - geHsH in the high field region. A 
dtailed analysis of the ESR shift in the whole range of H will 
be given in a subsequent publication.''-' 
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